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Abstract 

The freeze out of the expanding systems, created in relativistic heavy ion 
collisions, is discussed. We combine Bjorken scenario with earlier developed 

freeze out equations into a unified model. The important feature of the 
proposed model is that physical freeze out is completely finished in a finite 

time, which can be varied from (freeze out hypersurface) to oo. The 
dependence of the post freeze out distribution function on this freeze out 
time will be studied. As an example model is completely solved and 

analyzed for the gas of pions. 

In the ultrarelativistic heavy ion collisions at RHIC the total number of 
the produced particles exceeds several thousands, therefore one can expect 
that the produced system behaves as a "matter" and generates collective 
effects. Indeed strong collective flow patterns have been measured at RHIC, 
which suggests that the hydrodynamical models are well justified during the 
intermediate stages of the reaction: from the time when local equilibrium 
is reached until the freeze out (FO), when the hydrodynamical description 
breaks down. 

In the recent works [1] it has been shown that the basic assumptions of 
the Boltzman Transport Equation (BTE) are not satisfied during FO and 
therefore its description has to be based on Modified BTE, suggested in [Tj. 
When the characteristic length scale, describing the change of the distribu- 
tion function, becomes smaller than mean free path (this always happens at 
late stages of the FO), then the basic assumptions of the BTE get violated, 
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and the expression for the collision integral has to be modified to follow the 
trajectories of colliding particles, what makes calculations much more compli- 
cated. In fact, in cascade models one follows the trajectories of the colliding 
particles, therefore what is effectively solved is not BTE, but Modified BTE. 
Once the necessity of the Modified BTE and at the same time the difficulty 
of its direct solution were realized, the simplified kinetic FO models became 
even more important for the understanding of the principal features of these 
phenomenon. 

In simulations FO is usually described in two extreme ways: either FO on a 
hypersurface with zero thickness, or FO described by volume emission model 
or hadron cascade, which require an infinite time and space for a complete 
FO. However recently a new type of FO models has been developed in Refs. 
0, 13], in which freeze out is gradual and completely finished in a finite layer. 
Thickness of this FO layer can be varied from (freeze out hypersurface) 
to oo, and the dependence of the post freeze out distribution function on 
this freeze out time will be studied, making, thus, a bridge between the two 
extreme FO schemes discussed above. 

The model in Refs. [21 [3] is a further development of simplified FO models 
studied in Refs. jU, [3] , but in all these papers the expansion of the system was 
neglected. Thus, there is an important question to be studied: whether the 
important freeze out features, pointed out in these works, are not smeared 
out by the expansion of the system? 

In this paper we present a simple kinetic FO model, which describes the 
freeze out of particles from a Bjorken expanding fireball. Thus, this is a 
more physical extension of the oversimplified FO models without expansion 
|H 13, 121 13]- Taking the basic ingredients of the FO simulations from Refs. 
0, [3] we can make gradual freeze out to be completely finished in a layer, i.e. 
in a domain restricted by two parallel hypersurfaces r = T\ and r = T\ + L, 
where r is the proper time, r = \/t 2 — x 2 P, Ej. Thus, the evolution of the 
fireball created in relativistic heavy ion collision goes as follows. 
Initial state, r = tq: eo, tiq 

Phase I, Pure Bjorken hydrodynamics, To < r < t\ The evolution of 
the energy density and baryon density is given by the following equations [S] : 

de/r = — (e + P) jr , dn/dr = —n/r, (1) 

where P is the pressure. Bjorken model describes ID directed process - only 
the proper time gradients are considered, changes in the other directions are 
neglected. This system can be easily solved: 

e(r) = e (r /r) 1+c ° , n(r) = n (r /r) , (2) 
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where P = c 2 e is the equation of state (EoS) in general form. In order to have 
have a finite volume fireball, we need to put some boarders on the system. 
Here we assume that our system, described by the Bjorken model, is situated 
in the spatial domain \rj\ < ijr or what is the same \z\ < zr(t) = Tsmhrjn 
(r) = | In (frf) is pseudorapidity). As the system expends the volume of the 
fireball increases as V(r) = 2A xy sinh t]rt , where A xy is the transverse area 
of the system. 

Phase II, Bjorken expansion and gradual FO, t\ < r < t\ + L 

In order to describe freeze out we introduce two components of the distri- 
bution function, /: the interacting, /*, and the free (frozen out), f* ones, 
(/ = P + / )■ Now we can generate the interacting and free parts of baryon 
flow and energy-momentum tensor based on the corresponding components 
of the distribution function, which are parameterized in terms of interacting 
and free energy and baryon densities correspondingly. 

According to Refs. P, [7] the evolution of such a system is govern by the 
following equations for interacting and free components: 



dr t tfo L + T\ — r dr r rpo L + T\ — r 



(3) 



de? e l L drJ n l L 

— = + , -3— = + • 4) 

dr r tfo L + t\ — t dr r Tpo L + t\ — t 

In all these equations on the r.h.s. there are two terms: the first one is due to 

expansion (same as in Bjroken eqs. (CD)) and the second one is due to freeze 

out, see Refs. P, [7] for more details. 

Solving eqs. ([3]) we obtain: 

e (r) = e - , n l r x 



t) V L J 7 x y \tJ \ L J 

(5) 

The difference with respect to the pure Bjorken solution (J2j) is in the last 
multiplier, and we see that the interacting component completely disappears 
when r reaches the outer edge of FO layer r = L + t\. 

With these last equations we have completely determined evolution of the 
interacting component P, [7]. Knowing e l (r) and EoS we can find temper- 
ature, Ti(r). Due to symmetry of the system uf(r) = u^(tq) = (1,0,0,0). 
Finally, /*(r) is a thermal distribution with given 7^(r), n*(r), wf (r). 

However, the most interesting for us is the free component, which is the 
source of the observables. Eqs. (Ill) give us the evolution of the e/ and n/, and 
one can easily check that these two equations are equivalent to the following 
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equation on the distribution function 

dff 



dr 



T T F o L + T\ — T 



(6) 



The measured post FO spectrum is given by the distribution function at the 
outer edge of the FO layer, i.e. by f* (L + n). 

Most of the observables will depend only on this momentum distribution. 
However the two particle correlations depend also on the space-time origins of 
the detected particles. Thus, in order to calculate two particle correlations we 
have to keep track on both the momentum and the space-time coordinates of 
the freeze out point of particles (this can be done, for example, in the source 
function formalism), i.e. we need to know full evolution of free component 
and correspondingly can get some restrictions on it from the data. 

t 1 = 5.0 fm/c, x FO =0.5 fm/c 
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Figure 1: Evolution of temperature of the interacting matter for different FO layers. T,(ro = 
0.05 fm) = 820 MeV, T FO = 180 MeV. "No Freeze Out" means that we used standard 
Bjorken hydrodynamics even in phase II. 

To illustrate the FO process more quantitatively we show below the results 
for the ideal massive pion gas with Jiittner equilibrated distribution: 
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f(r J p} = 



(27T)* 



(7) 



where the degeneracy of pion is g = 3. 

Contrary to the illustrative example in [Bj here we do not neglect the pion 
mass. During FO the temperature of the interacting component decreases to 
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zero, so at late stages of the FO process this new calculation is better justified. 
We will see below that falls below m n quite soon, and so the Jiittner 
distribution is a good approximation of the proper Bose pion distribution [7]. 
For our system we have the following EoS: 
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m 2 T?K 2 (a) 



2tt< 



(8) 



where K n is Bessel function of the second kind, and a = m/TJ. 

Eqs. (J3j) and (JSj) result into the following equation for the evolution of the 
temperature of the interacting component: 

dTi Ti ^T?K 2 (a)+mT i K 1 (a) 



dr t 12T 2 K 2 {a) + bmT.K^a) + m 2 K {a) 



L 



ST 2 K 2 (a) + mT l K l (a 



r F0 \L + n- t) \2TfK 2 {a) + hmT.K^a) + m 2 K (a) ' v J 
Furthermore, we have used the following values of the parameters: t)r = 
4.38, A xy = kR\ u1 where Ra u = 7.685 fm is the Au radius, tq = 0.05 fm, 
T^tq) = 820 MeV , n = 5 fm, what leads to T^n) = T FO = 180 MeV, and 
Tpo = 0.5 fm. During the pure Bjorken case the evolution of the temperature 
is govern by eq. (jSj) without the second (freeze out) term on the r.h.s. 

In Fig. Q] we present the evolution of the temperature of the interacting 
matter, Tj(r), for different values of FO time L. 

x 1 = 5.0 fm/c, t fo =0.5 fm/c 
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Figure 2: Final post FO distribution for different FO layers as a function of the momentum 
in the FO direction, p = p x in our case (p y = p z = 0). The initial conditions are specified in 
the text. Dashed curve shows thermal distribution with temperature T = 160 MeV. 

As it was already shown in Ref. P, [3] , the final post FO particle distri- 
butions, shown on Fig. El are non-equilibrated distributions, which deviate 
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from thermal ones particularly in the low momentum region. By introducing 
and varying the thickness of the FO layer, L, we are strongly affecting the 
evolution of the interacting component, see Fig. but the final post FO 
distributions show strong universality: for the FO layers with a thickness 
of several tfo, the post FO distribution already looks very close to that for 
an infinitely long FO calculations, see Fig. E] left plot. Differences can be 
observed only for the very small momenta, as shown in Fig. [2] right plot. So, 
the inclusion of the expansion into our consideration does not smear out this 
very important feature of the gradual FO. 

Please note, that if one would look on our post FO distributions only in 
the medium momenta regions (with 0.5 GeV < \p\ < 2.5 GeV), then one 
could fit these spectra reasonable well with equilibrated distribution with 
temperature T = 160 MeV, see dashed line on Fig. El However, for low and 
high momenta such a fit would strongly disagree. 

Figs. [3] present the evolution of the total entropy, S(t), calculated based 
on the full distribution function, f(p) = f l (p) J r f^(p)' 



s(t) = / d s pf(r) 
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S(t) = s(r)V(r) 



(10) 



During pure Bjorken phase total entropy remains constant, as expected, but 
during phase II it constantly increases until FO is finished. 
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Figure 3: Evolution of the total entropy for different FO layers. The initial conditions are 
specified in the text. 



This is a very important conclusion of our model, which is stressing once 
again the importance to always check the non-decreasing entropy condition 



6 



P, El Ej, since long gradual freeze out may produce substantial amount of 
entropy, see Fig. [31 

This may have important consequence for QGP search. For many qual- 
itative estimations it was assumed that all the entropy is produced at the 
early stages of the reaction and that the expansion, hadronization and FO 
go adiabatically, and thus number of pions can serve as a rough measure of 
entropy. However, if a non-negligible part of the entropy, say 10% is produced 
during FO, then some estimations, for example of strangeness versus entropy 
(pion) production [TU] have to be reviewed. 
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